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INTRODUCTION 


The  mission  of  Che  GPS  control  segment  is  to  collect  measurements  from  remote 
unmanned  monitor  stations  at  a  central  location  -  the  master  control  station  (MCS)  - 
to  process  these  data  in  near  real  time,  and  to  upload  ephemeris  and  clock  data 
based  on  this  processing  to  each  of  the  GPS  satellites.  Figure  1  is  a  generalized 
graph  (timewise)  of  the  visibility  of  each  satellite  to  each  of  the  four  monitor 
stations.  Whenever  a  satellite  is  visible  to  one  or  more  monitor  stations,  data  is 
being  collected  and  transmitted  via  dedicated  phone  lines  to  the  MCS.  A  navigation 
message  can  be  produced  at  any  time  for  any  satellite  or  set  of  satellites  and  then 
uploaded  at  any  time  the  satellite(s)  are  visible  to  the  upload  station  (ULS).  The 
ULS  is  colocated  with  the  MCS.  In  practice,  daily  uploads  are  made  nominally  once 
a  day  at  a  time  to  provide  optimum  support  for  user  testing  at  the  Yuma  Proving 
Grounds.  Satellite  visibility  at  Yuma  is  about  the  same  as  at  Vandenberg  monitor 
station  (VMS),  so  a  study  of  figure  1  will  show  that  only  about  1.75  hours  of  four- 
satellite  visibility  is  available  at  Yuma  for  this  testing  each  day.  The  start  of 
this  period  is  the  rise  of  NAV  5.  The  desire  to  process  some  NAV  5  data  prior  to 
generating  the  upload  conflicts  with  the  need  to  provide  the  maximum  testing  time. 

A  compromise  is  made,  and  a  minimum  of  fifteen  minutes  of  NAV  5  data  is  processed 
prior  to  generating  the  navigation  message  for  this  satellite.  The  time  required 
to  perform  the  calculations  and  prepare  the  navigation  message  for  transmission  is 
nominally  3.5  minutes,  and  another  2.5  minutes  is  required  to  transmit  the  message. 
Yuma  ends  up  with  about  1.25  hours  of  test  time.  Clearly,  the  timeline  for  daily 
preparation  and  transmission  of  navigation  messages  is  tight.  Uploads,  during  the 
past  12  month  period  were  successfully  completed  according  to  this  schedule  97.257. 
of  the  time. 

Before  proceeding  through  the  total  sequence  of  operations  involved  in  the 
processing  of  data,  several  parts  of  the  total  process  will  be  discussed  separately 
to  avoid  lengthy  explanations  at  that  time. 


DATA  CPU  ACTION  SCHEDULE 

The  specific  collection  of  data  at  each  monitor  station  (MS)  is  controlled  via 
routine  messages  sent  to  each  MS  front  the  MCS.  While  it  is  possible  to  manually 
intervene  with  specific  instructions,  generally  routine  messages  assure  that  all 
available  data  is  collected  when  one  to  four  satellites  are  visible  at  a  MS.  When 
more  than  four  are  visible,  an  algorithm  that  resides  in  the  MS  computer  is  invoked 
to  time-share  the  MS  receiver's  (X-Set)  four  channels  with  the  visible  satellites. 

The  period  between  channel  switches  is  selectable,  and  currently  set  at  5  minutes, 
which  assures  that  with  up  to  six  satellites,  data  will  be  collected  on  all  satellites 
within  any  15  minute  interval. 

REFERENCE  EPHEMERIS 


The  Kalman  filter  employed  at  the  MCS  to  process  the  measurement  data  is 
linearized  about  a  reference  ephemeris  produced  by  the  Naval  Surface  Weapons  Center 
(NSWC)  for  each  satellite.  The  reference  ephemerides  are  produced  by  integrating 
40  day  trajectories  from  a  set  of  initial  conditions  based  on  a  least  squares  fit 
to  smoothed  GPS  measurement  data  (nominally  one  week's  worth)  transmitted  to  NSWC 
daily  by  the  MCS. 

tc 

NSWC  mails  a  new  set  of  reference  ephemerides  by  the  MCS  every  two  weeks.  The 
Kalman  filter  at  the  MCS  produces  residual  state  estimates  to  the  reference  values 
for  orbital  elements,  solar  pressure,  and  polar  wander. 

While  the  residual  estimates  produced  by  the  MCS  are  not  critically  dependent 
on  the  accuracy  of  the  predicted  reference  trajectories,  providing  linearity  is 
preserved,  they  are  the  fundamental  measure  of  the  fidelity  of  the  force  models 
being  used.  Considerable  effort  has  been  expended  to  isolate  sources  of  error 
in  this  respect.  A  major  improvement  was  realized  by  modification  of  the  solar 
pressure  model  as  a  result  of  work  carried  out  primarily  by  Mr.  William  Feess  of 


reference  trajectories  typically  run  off  in  the  along  track  direction  (which  is  by 

/ 

far  the  worst  case)  by  5-20  meters  per  day.  This  is  an  order  of  magnitude  improve¬ 
ment  over  what  was  originally  experienced.  Improved  data  from  the  monitor  stations, 
improved  processing  techniques  at  NSWC,  and  improved  solar  pressure  modeling  are 
each  responsible  in  part  for  this  improvement. 

SMOOTHING 

Given  a  range  and  delta  range  measurement  from  up  to  four  satellites  at  each 
of  four  monitor  stations,  from  a  practical  point  of  view,  it  is  necessary  to  compact 
these  total  data  for  input  to  the  Kalman  filter.  This  is  accomplished  by  polynomial 
fitting  the  range  and  delta  range  data  for  each  MS/SV  pair  over  a  fifteen  minute 
period  and  obtaining  a  single  pair  of  measurement  values  (range  and  delta  range)  for 
each  MS/SV  for  which  there  is  data.  The  Kalman  filter  is  cycled  once  each  fifteen 
minutes  providing  a  new  set  of  estimates  based  on  these  smoothed  data  pairs.  As 
will  be  clear  later  on,  the  smoothing  is  actually  done  to  corrected  measurement 
residuals.  For  the  purpose  of  obtaining  the  fit  for  each  15  minute  interval,  the 
range  and  delta  range  measurement  noise  is  assumed  uncorrelated,  and  weights  of 
.04  and  1  are  applied  for  the  range  and  delta  range  respectively.  Figures  2  and  3 
are  a  typical  plots  of  range  and  delta  range  residuals  respectively  to  the  fit  over 
a  15  minute  interval.  The  fit  is  linear  although  pro-vision  is  made  for  use  of  up 
to  a  fifth  order  polynomial.  The  fit  is  an  iteritive  least  squares  with  wild-point 
editing  incorporated.  Additionally  the  entire  15  minute  smoothed  value  is  rejected 
if  the  final  sample  standard  deviation  exceeds  a  selectable  value.  Nominally  this 
is  a  2ff#edit  applied  globally  without  regard  to  SV  or  MS  involved. 


KALMAN  FILTER  -  GENERAL 


The  MCS  Kalman  filter  employs  the  Carlson  square  root  algorithm  and  estimates 
6  residual  orbital  element  epoch  states,  two  epoch  solar  pressure  states  and  three 
epoch  clock  states  for  each  satellite.  In  addition  it  estimates  two  clock  states 
for  each  non-master  monitor  station,  one  tropospheric  scaling  state  for  each  monitor 
station  and  three  polar  wander  states.  The  filter  is  linearized  about  a  reference 
ephemeris.  With  the  range  and  delta  range  predicted  from  the  reference  subtracted 
from  the  measurements,  the  result  is  a  measurement  residual  (^R)  which  is  related 
to  the  residual  states  to  be  estimated  (S  X«)  by  the  linear  equation: 
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Ihe  partials  o  that  relate  to  the  ephemeris,  solar  pressure,  and 

polar  wander  states  are  provided  as  part  of  the  reference  ephemeris  by  NSWC.  The 
non-satellite  related  states  are  estimated  redundantly  in  each  active  partition. 
There  can  be  up  to  four  partitions  with  one  to  six  satellites  assigned  to  each.  A 
satellite  cannot  be  assigned  to  more  than  one  partition.  Each  partition  is  a 
totally  independent  filter.  As  a  consequence  of  a  data  transformation  performed  to 
decorrelate^  the  measurement  noises  (range  and  delta  range),  each  MS/SV  measurement 
pair  can  be  processed  sequentially  within  each  15  minute  Kalman  cycle.  Given  6 
satellites  and  4  monitor  stations,  a  maximun  of  24  smoothed  measurement  pairs  would 
be  processed  each  Kalman  cycle  if  all  satellites  were  simultaneously  visible  to 
all  monitor  stations. 


GPS  TIME 


User’s  of  the  GPS  system  determine  the  relative  range  to  each  of  4  satellites 


by  determining  the  time  offset  between  locally  generated  PRN  code  sequences  and 
identical  code  sequences  generated  and  transmitted  by  the  satellites.  These  code 
sequences  are  in  fact  the  output  of  the  respective  clocks.  The  navigation  message 
contains  polynomial  coefficients  that  predict  each  satellite's  clock  offset  relative 
to  a  common  time  scale.  This  time  scale  is  GPS  time. 

The  MCS  Kalman  filter  estimates  three  clock  states  for  each  satellite,  namely, 
time  offset,  and  its  first  and  second  derivitives;  or  equivalently  phase,  frequency 
and  frequency  drift.  These  estimates  are  with  respect  to  GPS  time.  The  filter  also 
estimates  2  states  (phase,  frequency)  for  three  of  the  four  monitor  station  clocks. 
The  fourth  monitor  station  contains  the  "master  clock".  No  estimates  are  made  for 
this  clock  as  it  is  the  time  base  for  the  GPS  system.  Any  monitor  station  can  be 
designated  as  the  "master".  Originally  the  Vandenberg  monitor  station  was  designated 
as  master,  and  GPS  time  was  the  same  as  the  VMS  station  clock.  Practical  considera¬ 
tions  have  required  master  clock  switches  from  time  to  time,  and  when  this  is  done, 
step  changes  in  GPS  are  avoided  by  bookkeeping  the  time  difference  between  old  and 
new  master  clock.  GPS  time  is  then  given  by  the  master  clock  adjusted  by  a  fixed 
time  bias  maintained  in  the  GPS  data  base.  This  bias  is  fixed  until  another  master 
clock  switch  is  made.  No  provision  is  made  to  prevent  a  step  change  in  the  rate  of 
GPS  time  when  a  master  clock  switch  is  made.  The  magr^tude  of  this  rate  change  of  •<_ 
course  depends  on  the  relative  frequency  offset  between  the  cesium  oscillators  in 
the  old  and  new  master  stations.  Typically  this  is  in  the  order  of  several  parts 
in  1013. 

PROCESSING  SEQUENCE 

Refer  to  figure  4  vhich  is  a  flow  chart  of  the  sequence  of  processing  steps 
described  below.  All  processing  as  described  takes  place  within  the  MCS. 


Raw  data  arrives  from  each  monitor  station  to  which  satellites  are  visible  as 
polled  by  the  MCS.  The  data  consists  of  pseudorange  and  delta  pseudorange  measure- 


merits  taken  each  6  seconds  and  time-tagged  by  the  monitor  station  clocks.  Differential 
range  (I4-L2)  to  be  used  for  ionospheric  correction  and  various  flags  indicating 
receiver  status  are  also  transmitted.  A  raw  data  edit  based  on  the  receiver  status 
flags  eliminates  any  meaningless  or  suspect  data  at  the  outset.  Prior  to  this  edit, 
a  magnetic  tape  (RAWRAN)  is  written,  from  which  if  necessary,  reprocessing  can  be 
initiated  in  non- real  time.  The  time  tags  are  corrected  to  GPS  time  of  signal  origi¬ 
nation  at  the  satellite.  Transit  time  and  monitor  station  clock  offset  is  determined 
by  dividing  the  pseudorange  measurement  by  the  velocity  of  light,  and  an  additional 
adjustment  for  satellite  clock  offset  is  made  based  on  the  current  Kalman  estimate 
of  this  state. 

A  reference  range  and  delta  range  pair  is  computed  by  interpolating  15  minute 
data  points  provided  in  the  reference  ephemeris  for  the  corrected  GPS  time  tag  of 
each  measurement.  These  reference  values  are  subtracted  from  the  measurements  to 
generate  residual  measurements.  These  residual  measurements  are  the  basis  for 
estimating  residual  states  in  the  Kalman  filter.  These  residuals  are  then  corrected 
for  deterministic  sources  of  measurement  error  as  listed  in  figure  4. 

At  this  point,  corrected  measurement  residuals  for  each  6  second  range  and  delta 
range  measurement  that  passed  the  raw  data  edit  for  all  active  MS/SV  pairs  are 
accumulated  in  an  ordered  file  for  a  given  15  minute  interval.  When  all  available 
data  are  collected  for  a  given  interval,  the  smoothing  algorithm  is  invoked  to 
create  a  single  pair  (range  and  delta  range)  of  smoothed  measurements  for  each  MS/SV 
pair  for  that  15  minute  interval.  The  smoothing  algorithm  is  a  weighed  least  squares 
polynomial  (linear)  fit  to  the  combined  range  and  delta  range  data,  with  the  delta 
range  weighted  25  times  the  range  data.  Individual  wild-point  editing  is  performed 
with  iteration,  and  the  final  values  are  edited  based  on  the  sample  variance  {2c*  ). 
Each  pair  of  smoothed  values,  is  edited  again  by  comparison  with  an  expected  value 
and  its  uncertainty  based  on  the  last  set  of  Kalman  estimates  and  their  covariances. 

Not  shown  in  the  flow  chart,  is  a  transformation  applied  to  the  smoothed 


i  yank m  in 


[Jl- 


measurement  residual  pairs  to  decorrelate  the  range  and  delta  range.  (In  the  smoothing 
algorithm,  the  smoothed  value  selected  lies  midway  between  the  earliest  and  lastest 
data  point.  If  this  is  not  exactly  midway  in  the  15  minute  interval,  some  correla¬ 
tion  is  introduced.)  Decorrelation  at  this  point  permits  sequential  processing  of 
the  smoothed  measurement  residuals  in  the  Kalman  filter. 

The  Kalman  filter  updates  the  estimates  of  its  residual  states  once  each  15 
minutes,  and  at  any  time,  a  prediction  of  each  satellite's  position  and  clock  off¬ 
set  can  be  generated  fran  these  residual  estimates  together  with  the  reference 
values.  Refer  to  figure  5  for  the  remainder  of  this  processing  sequence. 

Given  Kalman  residual  state  estimates,  the  satellite  position  can  be  predicted 
in  cartesian  earth-fixed  coordinates  by  the  equation  shown  in  figure  5.  To  generate 
a  navigation  message  say  for  26  hours  into  the  future,  this  equation  is  employed  to 
generate  26  ninety  minute  sets  of  cartesian  satellite  coordinates,  and  a  least 
squares  fit  is  made  for  a  set  of  14  Keplerian-type  elements  to  each  ninety  minutes. 
There  is  a  half  hour  overlap  of  each  set,  so  that  a  new  set  is  to  be  used  for  each 
hour.  The  writer  has  employed  the  standard  user  algorithm  in  a  hand-held  Tl-59 
calculator  to  compute  satellite  cartesian  coordinates  from  these  elements,  and 
obtained  agreement  within  a  few  centimeters  with  the  cartesian  coordinates  which 
formed  the  basis  for  creating  the  elements. 

The  satellite  clock  predictions  are  made  by  propagating  the  individual  epoch 
clock  state  vector  (  &Te  t  &/T0  L Te)  by  the  equation  shown  in  figure  5,  and  calcul¬ 
ating  a  set  of  coefficients  (  &0  ,  Q_ «  ,  0,  )  and  a  reference  time  A.  ,  for  each 

hour  (ninety  minutes  with  overlap).  The  clock  coefficients  are  each  perturbed  by 
a  function  vAdch  describes  the  effect  of  general  relativity  so  that  users  may 
ignore  this  effect  in  their  measurements. 


SYS'TM  PERFORMANCE  MEASURES 


STATE  ESTIMATES 


Plots  are  made  of  selected  residual  state  estimates  on  a  routine  basis.  For 
example  figure  6  is  a  30  hour  plot  of  the  Kalman  estimate  of  current  time  frequency 
for  NAV  3  -  a  rubidium  Figure  7  is  the  same  for  NAV  5  -  a  cesium.  The  straight 
line  portions  of  the  curve  are  during  periods  of  non-visibility  of  the  satellite, 
and  represent  the  propagated  value.  Note  the  slope  of  this  straight  line  portion 
for  the  rubidiun  (NAV  3)  indicating  a  nonzero  estimate  for  the  frequency  drift 
state.  Excessive  variation  in  estimates  of  clock  states  can  be  an  indication  of 
ephemeris-clock  state  estimate  aliasing.  When  this  occurs ,  reinitialisation  of 
the  filter  covariance  matrix  is  used  to  correct  the  problem. 

USER  RANGE  ERROR 


As  an  overall  measure  of  system  performance,  a  quantity  called  user  range 
error  (URE)  is  computed  and  plotted  daily.  The  URE  is  a  measure  of  the  difference 
between  ephemeris  and  clock  predictions  currently  active  in  the  navigation  message, 
with  tho'-'  same  variables  according  to  current  Kalman  estimates.  The  assumption  is 
that  current  Kalman  estimates  most  accurately  represent  the  truth,  and  the  difference 
is  a  measure  of  the  error  in  the  navigation  message.  The  ephemeris  part  of  the  URE  is 
that  portion  of  the  total  satellite  position  error  projected  along  the  line-of-sight 
to  a  selected  user  location  (nominally  Yuma  Proving  Grounds). 

Figure  6  is  a  "good”  URE  and  figure  7  is  a  "poor"  URE.  Almost  all  of  the  time, 
URE's  range  between  these  limits  with  the  greater  percentage  looking  like  the  former. 
The  "total  URE"  is  the  net  resulting  range  error  due  to  both  the  ephemeris  and  clock 
errors  in  the  navigation  message.  The  run-off  in  the  ephemeris  component  as  seen  in 


figure  7  is  for  the  most  part  the  geometric  effect  of  a  more  or  less  fixed  along- 
track  ephemeris  error.  Near  the  end  of  the  plot,  the  satellite  is  near  setting. 

YUMA  REPORTS 

Daily  measures  of  system  performance  are  made  by  the  test  facility  at  the  YPG. 
They  ccmpute  relative  satellite  clock  synchronization  errors  from  range  and  delta 
range  measurements  at  a  known  location  and  also  measure  static  navigation  errors. 
Generally  these  measures  correlate  qualitatively  on  a  day  to  day  basis  with  MCS 
measures . 

ANALYTICAL  DESCRIPTION  OF  MCS  KALMAN  FILTER 

The  basic  model  for  the  GPS  Kalman  Filter  is 

=  <|  6  K+i ,  to  X.(K')  +  wdk+rt 
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These  equations  apply  to  estimates  for  vehicles  within  a  given  partition, 
and  of  course  the  number  of  states  is  a  function  of  the  number  of  vehicles  repre¬ 
sented  within  a  partition.  Each  vehicle  is  represented  in  a  twelve  dimensional 
vector  space,  the  first  six  states  representing  Keplerian  orbital  elements,  the 
next  three  solar  pressure  (two  used)  states,  and  the  last  three  clock  terms.  There 
is  also  a  thirteen  dimensional  monitor-station  state  space  to  represent  the  four 
monitor-stations  currently  used  by  GPS  for  data  collection.  There  are  two  clock 
states  for  each  monitor  station  except  a  distinguished  station,  designated  master 


a-posteriori  covariance  from  the  fitting  process.  \ 

The  measurement  update  for  the  system  is  K  (  2  '  W  *  )  where  Z  is  vector  of 

the  form  (  2  it  ,Z2>,- •  •  •  2  r2.,2*i” (i.e. ,  the  processing  is  done 

through  vehicles  for  a  fixed  monitor  station),  each  2;^-  ( ,ATt^  ,  K.  being  the 
gain  computed  as  K  =  (HRit  +R)-^,  with  P  the  covariance  updated  through  the  time 

update,  and  H  being  the  measurement  matrix  "H  =  (  Wn  •••WSI  Un, . VLj) ,  each 

H  i  ^  f  •  j  with  the  measurement  matrix  for  the  ith  vehicle 

and  jth  monitor  station.  (This  is  the  standard  algorithm  filtering.  Actually,  the 
GPS  system  uses  the  Carlson  square  root  algorithm,  but  for  our  purposes  here  this 
is  unessential). 

Each  Ht'i  is  a  function  of  time  tags  based  on  the  data  smoothing  during  the 
measurement  interval.  If  "t  ka  ■ .  ,t.E*  ‘  >  and t- •  represent  the  mid- point,  end- 
point  and  start  point  of  the  smoothing  interval,  then  Ht  j  rH4'  ^  (i  J  =  -H(tsY 

Furthermore,  each  of  the  measurement  matrices  can  be  factored  as  t  pj  where 

-iw  1 

relates  only  to  vehicle  states  and  H  j  only  to  the  jth  monitor  station  states. 

Because  of  the  assumption  that  all  measurements  between  vehicle  and  monitor 

j 

station  are  independent,  the  measurement  update  can  be  processed  sequentially  by 
satellite-monitor  station  pair.  This  goes  as  follows: 

If  the  update  were  done  as  a  full  matrix  process  (ignoring  the  square-root 
algorithm)  the  update  actually  processed  is: 

X*  +  K 

where  r  P  WT  (  W  ?  HT  *  2  ) 


Thus,  each  measurement  update  only  involves  the  matrix  .  .  Using  this 

*-  *■  i 

estimate  of  the  measurement  noise  for  the  (  c  j  )  data  pair,  a  decorrelation  of 

range  and  delta-range  is  accomplished  by  transforming  to  variables 2**  .Housing 

/w  j/2  _  ~  -  Vz  -  -yz  * 

:  T’-  S-tj  f  M  *  H  ,  where  £  •  •  is  che  lower-diagonal 


Cholesky  square  root  of  S  t ’  •  If  ^ 


r« 


then  Y 


^  "  \f«<  J  ‘  -zL_ 

L  J  (5*6-p*){ 

The  decorrelated  variables  have  a  measurement  noise  covariance  equal  to  2  x  2 
identity  matrix,  which  allows  the  range  and  delta-range  for  a  given  vehicle-monitor 
station  pair  (  \  )  to  be  processed  sequentially.  (For  more  information  regarding 

the  technique  of  decorrelating,  as  well  as  details  regarding  various  square-root 
algorithms  and  sequential  processing,  see  Factorization  Methods  for  Discrete 
Sequential  Estimation,  Gerald  J.  Bierman,  Academic  Press,  1977.) 

Assuming  the  decorrelations  have  been  carried  out,  the  sequential  processing 
proceeds  by  cycling  through  the  following  set  of  equations  for  each  vehicle-monitor 
station  data  pair:  KllJ+ri*  Pi‘-»  Ul/(V\L  Pi~%  \n  i  +  l) 

n.i  -•  Czi-i’u 

Pi-. 

where  PQ  and  Xq  are  the  a-priori  covariances  and  states  after  time  update ,  "2  i  rep¬ 
resent  measured  range  and  measured  delta-range  (processed  one  after  the  other)  and 
the  final  value  update  after  processing  all  available  data  is  the  value  used  in  the 
next  time-update.  Note  that  each  gain  K  'L  has  a  scalar  in  its  denominator.  Also 
note  that  this  update  gives  a  new  estimate  of  the  total  state  due  to  the  ( 
range  and  delta-range.  Processing,  as  noted  before,  cycles  through  vehicles  per 
monitor  station. 

Another  way  to  view  the  processing  is  by  examining  the  measurement  update  process 
globally.  Using  the  partitioning  of  the  state  space  outlined  above,  PHC  can  be  cal¬ 
culated  as  a  matrix  with  12  x  2  matrix  coefficients,  further  subdivided  by  monitor 
stations,  so  that  PHC  =  [PH^,  PH^,  PH§,  PHE],  with  each  PH£  of  the  form 

(Ph  Hs<  PimHmHPiiHs  ♦  P i ►a H •  •  '(Pi's Ha  4 Piw H 
(Pl2  VW  ♦  PzHvita)  Hs  *  •  *(Pfc*  TU  ♦  PaM  Hi*) 

(^5«  Hi  4PsM  HwJ  (PszHs  *  Psm  H  *  '  *(Ps5  iU  +  P5*5l*) 

Ms  4  H  tPlAIAHl*)'  *  *(PnsHS4  Pia*a  ! 


station,  for  which  no  clock  states  are  estimated  (its  clock  becomes  the  system 
clock).  There  is  also  a  state  for  a  tropospheric  correction  at  each  monitor  sta¬ 
tion  and  three  polar  wander  states  which  are  cannon  to  all  the  monitor  stations. 


This  means  that  a  partition  represents  a  state-space  of  dimension  12S  +  13, 
vhere  S  is  the  number  of  vehicles  in  a  partition.  In  the  GPS  system,  1  £  5  £  4  . 

The  estimated  states  are  residual  corrections.  The  ephemeris  portion  is 
linearized  about  a  nominal  reference  trajectory,  clock  states  are  relative  to  GPS 
time,  and  all  states  are  referenced  to  an  epoch  time  (which  is  different  for  clock 
and  ephemeris  states). 

If  V  represents  the  state  space  for  a  single  partition,  then  is  decom¬ 
posed  into  a  direct  sum  of  subspaces  ,  with  one  sub-space  for  each 

vehicle  and  a  single  subspace  for  the  monitor  station  states. 

The  state  transition  matrix  <J>  is  blocked  to  acccnmodate  the  above  scheme  as 


o 


L  °  in  I 

Hlth  <$\)  -$,V,  ■  +  4m Vm  ,  each  4  i  leaving  V; 

invariant.  Because  of  this,  the  time  update  is  accomplished  independently  for  each 

yJl  .  However,  since  the  filter  is  an  epoch  state  filter,  the  process  noise  is 
^  •*  T 

calculated  via  P  -  p  Per  where  P  is  the  a-priori  covariance  matrix 


for  the  filtering  error  and  X  =  X,  X  being  the  a-priori  state  estimate. 

The  measurement  update  is  processed  for  range  and  delta-range  measurements 
between  a  given  vehicle  and  a  given  monitor  station.  The  two  measurements  are 
assumed  independent  and  are  collected  every  six  seconds  for  a  15  minute  interval. 
This  data  is  smoothed  through  a  linear  weighted  least-squares  estimate,  with  the 
weighting  matrix  a  fixed  diagonal  matrix  weighting  delta-range  five  times  as  heavily 
as  range.  The  output  of  this  process  is  a  pair  £***>  A^J^of  range  and  delta-range 
residuals  between  the  i**  vehicle  and  j*'1  MS,  together  with  *  t^e 


A  ^ 

with  each  of  Yi\  Us  •»  Pi*  W*  a  12  x  2  matrix  and  a  13  x  2  matrix. 

If  one  multiplies  the  first  row  (of  12  x  2  matrices)  times  the  measurement  vector 
(  (  t  ii  t  M“,« the  effect  will  be  the  total  update  of 


where  (  r*^  ^  }  k  )  is  the  measurement  residual ,  P  t'  •  is  the  submatrix  of  P 

•  ^  *  >  — 

representing  the  covariances  between  vehicle  {.5  states  and  vehicle  J  *  states,  P^K 
is  the  cross-correlation  matrix  for  the  i+v>  vehicle  and  the  monitor  station 
states.  From  this  calculation  one  can  determine  the  separate  effects  on  the  vehicle  t 
due  to  data  from  vehicle  \  and  monitor  station  K  through  the  covariance  matrix 
P.  By  stripping  off  the  appropriate  covariances  and  measurement  matrix  values  each 
of  these  computations  can  be  examined  individually.  The  term  P»j  ^5  relates  the  i** 
vehicle  states  to  the  i*1'  vehicle's  states,  while  relates  the  VC***  monitor 

station  states  to  the  vehicles  states. 


(The  last  "row"  of  13  x  2  matrices  produces  the  monitor  station  estimates  based 
on  all  the  data.) 


To  tie  this  in  with  the  sequential  processing,  notice  that  the  sequential 
update  operates  on  one  pair  (  )  at  a  time  (before  decorrelation,  when  the 

pairs  are  separated),  taking  a  column  of  dimension  two  in  width  and  multiplying  it 


by  the  appropriate  pair.  This  gives  the  effect  on  all  the  state  estimates  due  to  the 

^  /yj  ^  ^ 

appropriate  pair  (  Ct"  j  ^  ).  (For  example,  colunn  one  is  multiplied  by  (*“,* ,  aqN) , 

colunn  two  by  (^ji  ,  )  etc.)  Thus,  the  sequential  processing  can  be 


thought  of  as  colunn  processing, 
allow  sequential  processing  of 


is  to 


CONCLUSIONS 


Th*  Plus*  I  GPS  has  exhibited  excellent  performance  as  a  real-time  navigation 
system.  It  remains  to  be  seen  if  accuracies  equivalent  to  TUANS  IT  are  possible 
far  geodetic  applications.  To  a  large  extent  this  may  depend  an  the  development 
of  optimum  processing  strategies  with  multiple  passes  or  aultiple’ satellite 
data.  It  is  not  clear  at  this  time  ehether  precise  epheaerides  can  be  produced 
that  are  significantly  better  than  the  current  navigation  data  transmitted  from 
the  satellites. 
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